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Abstract 



We elaborate on a new technique for computing properties of nucleon-nucleon 
interactions in terms of an effective field theory derived from low energy NN 
scattering data. Details of how the expansion is carried out to higher orders 
are presented. Analytic formulae are given for the amplitude to subleading 
order in both the % and ^ - 3 £>! channels. 
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I. INTRODUCTION 



Effective field theory appears to be an ideal tool for the study of low energy nuclear 
physics, as the nucleon energies are typically well below the complex spectrum of hadrons 
that exist with masses greater than about lGeV. An example of the successful application 
of effective field theory to low energy hadronic physics is chiral perturbation theory, which 
exploits the fact that the lightest pseudoscalar mesons are approximate Goldstone bosons 
(for a recent review, see Even though an analytic description of pions in terms of quarks 
and gluons is impossible, our ignorance can be parametrized in an effective theory such that 
a perturbative calculation of pion interactions involving only a few parameters agrees well 
with experiment. Central to the utility of chiral perturbation theory for mesons is that there 
is a clear power counting scheme, so that one can include all effects to a given order, and 
estimate the size of errors incurred in the approximation. 

Recently we applied the same procedure to nucleon-nucleon interactions, outlining a 
method to consistently expand the NN interaction in powers of p and m n , where p is the 
momentum of each nucleon in the center of mass frame, and m n is the pion mass |§ . The 
analysis was inspired by Weinberg's proposal || that effective field theory could be profitably 
used in nuclear physics, as well as by subsequent work fl-[T7|. The original idea was to exploit 
the approximate chiral symmetry of the strong interactions, which gives rise to a hierarchy 
of length scales between the Compton wavelengths of the vector mesons and the pions. In an 
effective field theory, short distance nucleon-nucleon interactions are encoded in a derivative 
expansion of local operators. This is in contrast with the various models of extended nucleon- 
nucleon potentials with free parameters chosen to fit scattering data. These models can fit 
the nucleon-nucleon phase shifts to great accuracy, but suffer from several deficiencies not 
shared by the effective field theory approach: they are not useful for computing inelastic 
processes, they give no insight into three-nucleon forces, and they are numerically intensive 
to use in the iV-body problem. Furthermore, there is no systematic way of anticipating the 
errors one should expect when using these potentials. Another advantage of effective field 
theory is that it can easily incorporate chiral symmetry, and can be naturally extended to 
discuss systems with strange quarks, such as hypernuclei [JTH] and kaon condensation ]TP| , |2"U| |. 

However, the effective field theory analysis of the two-nucleon system is complicated by 
the existence of other length scales, in particular the S'-wave scattering lengths, which are 
many times longer than the pion Compton wavelength. The existence of large scattering 
lengths implies that the underlying physics at short distance is both nonperturbative and 
"finely tuned", in the sense that the interactions must be near a critical value. In Weinberg's 
original work || a power counting scheme was proposed that involved summing a particular 
infinite class of Feynman graphs at each order in the expansion. However, it was shown in 
refs. [0,|6| that graphs at a given order in that expansion require counterterms corresponding 
to operators treated as being higher order, and that therefore the expansion proposed in |3j 
is inconsistent. This was also demonstrated in models where a perturbative matching could 
be performed jl2| |. 

Ref. presented a different expansion that cures this problem, and applied it to NN 
scattering in the ^So and 3 S\ — 3 D± channels. In this paper we begin by giving further 
details about the expansion, analyzing in detail a model of heavy bosons interacting via 
meson exchange. This simplified theory serves to explain the special treatment accorded 
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systems with large scattering lengths, and explains the virtues of the PDS subtraction scheme 
and renormalization group analysis introduced in ref. 0. We then explain how the power 
counting is extended to the theory with pions and give the explicit analytic formulae for the 
spin-triplet scattering amplitudes to subleading order. 



II. EFFECTIVE FIELD THEORY FOR NONRELATIVISTIC SCATTERING: A 

TOY EXAMPLE 

In this section we present a toy model of heavy spinless "nucleons" N interacting via a 
Yukawa interaction characterized by a scale A. We then construct the effective field theory 
describing scattering at momenta p C A, consisting entirely of contact interactions in a 
derivative expansion. Since these local operators are singular, this formulation of the low 
energy theory necessarily introduces divergences that must be dealt with by the conventional 
regularization and renormalization procedures, so that the final result is independent of a 
momentum cutoff. We show how to organize the Feynman graphs in the effective theory in 
a consistent power counting scheme so that the scattering amplitude can be expanded in 
powers of p/A. Since the sizes of all the coupling constants in the effective theory depend 
on the subtraction scheme used to render diagrams finite, the development of the power 
counting scheme is intimately related to the renormalization procedure used. We show why 
the PDS subtraction scheme introduced in is particularly well suited for this problem, 
and we show how a renormalization group analysis is useful in the case of systems with large 
scattering length, such as those seen in the realistic problem of nucleon-nucleon scattering. 

It should be no surprise that the effective field theory expansion for the toy system is 
simply related to the conventional effective range expansion, and so the machinery of quan- 
tum field theory may appear to be heavy handed and superfluous. Nevertheless, the field 
theoretic language that we develop in this section is readily extended to the realistic problem 
of interest: nucleons interacting via both short range interactions and long range pion ex- 
change. In the realistic problem, effective field theory is not equivalent to an effective range 
expansion, and is the only framework that can consistently incorporate chiral symmetry and 
relativistic effects without resorting to phenomenological models. 

We assume that the spinless bosons iV are nonrelativistic with mass M, carry a conserved 
charge ("baryon number"), and interact via the exchange of a meson <fi with mass A and 
coupling g. At tree level, meson exchange gives rise to the Yukawa interaction 

a 2 e~ Ar 

V{r) = — , (2.1) 
47r r 

and the Schrodinger equation for this system may be written as 



„ e p 
x x A 2 



* = , (2.2) 



x = Af, 71= 9 -^-, p 2 = ME . (2.3) 
Arc A 

Note that p is the magnitude of the momentum carried by each N particle in the center of 
mass frame. Evidently there are two options for a perturbative solution for the S'-matrix 



3 



for this system. The first is an expansion in powers of rj, the familiar Born expansion. An 
alternative is to expand in powers of p/A, which is the expansion parameter used in effective 
field theory. An important feature of the low energy expansion is that it can provide accurate 
results in terms of a few phenomenological parameters even for nonperturbative rj. This is 
the regime we are interested in, and so we will assume throughout that rj ~ 1. 

The quantity that is natural to calculate in a field theory is the sum of Feynman graphs, 
which gives the amplitude iA, related to the S'-matrix by 

S = i + iMlA. (2.4) 
For S'-wave scattering, A is related to the phase shift 5 through the relation 

A=^ \ . (2.5) 

M p cot 5 - ip 

From quantum mechanics it is well known that it is not A, but rather the quantity p cot S, 
which has a nice momentum expansion for p <C A (the effective range expansion): 

1 1 00 /» 2 \ n+1 

where a is the scattering length, and r is the effective range. So long as rj ~ 1, the case 
we will be interested in, the coefficients r n are generally 0(1/ A) for all n, but a can take 
on any value, diverging as i] approaches one of the critical couplings r]k for which there is a 
boundstate at threshold. (The lowest critical coupling is found numerically to be r/i = 1.7.) 
Therefore the radius of convergence of a momentum expansion of A depends on the size 
of the scattering length a. In the next section the situation where the scattering length is 
of natural size \a\ ~ 1/A is considered, while in the subsequent section we discuss the case 
\a\ ^> 1/A, which is relevant for realistic NN scattering. 



A. The momentum expansion for a scattering length of natural size 



In the regime \a\ ~ 1/A and \r n \ ~ 1/A, A has a simple momentum expansion in terms 
of the low energy scattering data, 



A = 



Ana 



l-iap+ (ar /2 - a 2 )p 2 + 0(p i / 'A J ) 



(2.7) 



which converges up to momenta p ~ A. It is this expansion that we wish to reproduce in 
an effective field theory. 

The effective field theory of N particles interacting through contact interactions has the 
following Lagrangian: 

C = A>t (idt + V 2 /2M) N + {fi/2) 4 - D [c (A>tA>) 2 + cP(N^VN) 2 + C? ] \iV N)] 2 + .. 

The sum of Feynman diagrams computed in this theory gives us the amplitude A. As we 
will be using dimensional regularization for the loop integrals in this theory, the spacetime 



(2.8) 
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FIG. 1. The bubble chain arising from local operators. 



dimension is given by D. Dimensional regularization is the preferred regularization scheme as 
it preserves gauge symmetry and chiral symmetry, as well as Galilean invariance (or Lorentz 
invariance, for relativistic systems). The ellipses indicates higher derivative operators, and 
(fJ>/2) is an arbitrary mass scale introduced to allow the couplings C 2n multiplying operators 
containing V 2n to have the same dimension for any D. We focus on the S-wave channel 
(generalization to higher partial waves is straightforward), and assume that M is very large 
so that relativistic effects can be ignored. The tree level S'-wave amplitude is 

oo 

iA tiee = -*(/i/2) 4 " D £ C 2 Mp 2n , (2.9) 

n=0 

where the coefficients C 2 n{^) are various linear combinations of the couplings in the La- 
grangian eq. (|2.8f) contributing to S'-wave scattering. 

The loop integrals one encounters in diagrams shown in Fig. [I] are of the form 



(2tt) d ^ \E/2 + q - q 2 /2M + ie J \E/2 - q - q 2 /2M + ie, 



(27r)(°- 1 ) m \E-q 2 /M 



ie 



\A-D 



= -M(ME) n (-ME - i er-'^r (L_°) ^l 1)/2 • (2.10) 

In order to define the theory, one must specify a subtraction scheme; different subtraction 
schemes amount to a reshuffling between contributions from the vertices and contributions 
from the the UV part of the loop integration. For the case we are considering, \a\, |r n | ~ 1/A, 
it is convenient to use the minimal subtraction scheme (MS) which amounts to subtracting 
any 1/(D — 4) pole before taking the D — > 4 limit. The integral eq. (|2.10|) doesn't exhibit 
any such poles and so the result is simply 

I;?* = (MET (^) V-ME - ie = -i (^) . (2.11) 

A nice feature of this scheme is that the factors of q inside the loop get converted to factors of 
p, the external momentum. Therefore one can use the on-shell, tree level amplitude eq. ( f2.9|) 
as the internal vertices in loop diagrams and summing the bubble diagrams gives 

A = g ° 2np2n (2 12) 

l + i(Mp/Aix)Y.C 2n p^ ■ { ■ ! 



5 



Since there are no poles at D = 4 in the MS scheme the coefficients C2n are independent of 
the subtraction point \i. The power counting in the MS scheme is particularly simple: each 
vertex C^V 2 ™ counts as order p 2n , while each loop brings in a factor of p. The amplitude 
may be expanded in powers of p as 

oo 

A=J2 A n , A n ~0(p n ) (2.13) 

n=0 

where the A n each arise from graphs with L < n loops and can be equated to the low energy 
scattering data eq. (|2.7|) in order to fit the C*2 n couplings. In particular, Ao arises from the 
tree graph with Co at the vertex; At is given by the 1-loop diagram with two Co vertices; 
A2 is gets contributions from both the 2-loop diagram with three C vertices, as well as the 
tree diagram with one C2 vertex, and so forth. Thus the first three terms are 

•4o = -Co, A = <^, A 2 = Cl{^) -C 2 p 2 . (2.14) 



Comparing eqs. (|2.7| , |2.14j) we find for the first two couplings of the effective theory 



Co = ^, C 2 = C^. (2.15) 
In general, when the scattering length has natural size, 

°* ~ mik ■ < 2 ' 16) 

Note that the effective field theory calculation in this scheme is completely perturbative 
even though rj ~ 1 and there may be a boundstate well below threshold. The point is, 
that when there are no poles in A in the region \p\ < A, the amplitude is amenable to a 
Taylor expansion in p/A in that region, and with a suitable subtraction scheme this Taylor 
expansion can correspond to a perturbative sum of Feynman graphs. 



B. The momentum expansion for large scattering length 



Now consider the case \a\ 3> 1/A, \r n \ ~ 1/A, which is of relevance to realistic NN 
scattering. For a nonperturbative interaction (77 ~ 1) with a boundstate near threshold, the 
expansion of A in powers of p is of little practical value, as it breaks down for momenta 
p > l/\a\, far below A. In the above effective theory, this occurs because the couplings C2n 
are anomalously large, C2 n ~ 47ra n+1 /MA n . However, the problem is not with the effective 
field theory method, but rather with the subtraction scheme chosen. 

Instead of reproducing the expansion of the amplitude shown in eq. Q2.7Q , one needs to 
expand in powers of p/A while retaining ap to all orders: 



A 



Att 



1 



M (1/a + ip) 



1 + 



r /2 



(n/2A 2 

(1/a + ip)*' (l/a + ip) 2 *' ' (1/a + ip) 



p* + 



-p 4 + 



(2.17) 



Note that for p > l/\a\ the terms in this expansion scale as {p 1 ,p°,p 1 , . . .}. Therefore, the 
expansion in the effective theory should take the form 
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A=Y i A n , A n ~0(p n ) (2.18) 



instead of the expansion in eq. Q2.13|) . Again, the task is to compute the A n in the effective 



theory, equate to the appropriate expression in eq. ( |2.17| ), and thereby fix the C 2n coefficients. 
For example, 

As we have seen in the previous section, any single diagram computed in the effective 
theory is proportional to positive powers of p. The leading term A-i must therefore involve 
summing an infinite set of diagrams. It is easy to see that the leading term in eq. (|2.17 ) can 
be reproduced by the sum of bubble diagrams with C vertices H , which yields 

A-l = r "5?, n , (2.20) 



1 + 



in the MS scheme. Comparing this with eq. ( 2.19|) we find Co = Ana/M, as in the previous 



section. However, there is no expansion parameter that justifies this summation: each 
individual graph in the bubble sum goes as (4iia/M)(iap) L , where L is the number of loops. 
Therefore each graph in the bubble sum is bigger than the preceding one, for \ap\ > 1, while 
they sum up to something small. 

This is an unacceptable situation for an effective field theory; it is important to have an 
expansion parameter so that one can identify the order of any particular graph, and sum 
them up consistently. Without such an expansion parameter, one cannot determine the 
size of omitted contributions, and one can end up retaining certain graphs while dropping 
operators needed to renormalize those graphs. This results in a model-dependent description 
of the short distance physics, as opposed to a proper effective field theory calculation. 

Since the sizes of the contact interactions depend on the renormalization scheme one 
uses, the task becomes one of identifying the appropriate subtraction scheme that makes 
the power counting simple and manifest. The MS scheme fails on this point; however this is 
not a problem with dimensional regularization, as has been frequently suggested, but rather 
a problem with the minimal subtraction scheme itself. The momentum space subtraction at 
threshold used in ref. J| behaves similarly. 

An alternative regularization and renormalization scheme is to use a momentum cutoff 
equal to A. Then for large a, Co ~ (Att/MA), and each additional loop contributes a factor 
of Co(A + ip)M/4iT ~ (1 + ip/A). For A ^> p the factor of p/A is small relative to the 1, 
but neglecting it would fail to reproduce the desired result eq. ( j2.19| ). The problem is that 
there are significant cancellations between terms in this scheme. 

Evidently, since A-i scales as 1/p, the desired expansion would have each individual 
graph contributing to A-i scale as 1/p. As the tree level contribution is Co, we must have 
Co be of size cx 1/p, and each additional loop must be 0(1). This can be achieved by using 
dimensional regularization and the PDS (power divergence subtraction) scheme introduced 
in ref. 0. The PDS scheme involves subtracting from the dimensionally regulated loop 
integrals not only the 1/(D — 4) poles corresponding to log divergences, as in MS, but also 
poles in lower dimension which correspond to power law divergences at D = 4. The integral 
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I n in eq. ( |2.10|) has a pole in D = 3 dimensions which can be removed by adding to I n the 
counterterm 



M(ME) n n 



Att(D - 3) ' 

so that the subtracted integral in D = 4 dimensions is 



(2.21) 



tPDS 



In + SI n 



-{ME) 7 



M 

47T 



(2.22) 



An alternative subtraction scheme with similar power counting is to perform a momentum 
subtraction at p 2 = —fi 2 , as recently suggested in ref. pl|. The PDS scheme retains the 



nice feature of MS that powers of q inside the loop integration are effectively replaced by 
powers of the external momentum p. Note that at [i = 0, PDS and MS are the same. In 
this subtraction scheme 



A 



2n 



M(jJ + tp)/47Tj:C2nP 2n 



(2.23) 



The amplitude A is independent of the subtraction point fi and this fact determines the 
\i dependence of the coefficients, C*2 n - In the PDS scheme one finds that for \i 3> 1/H, the 
couplings C*2n(^) scale as 



C 2 „(//) 



47T 



MA n fi 



n+l 



(2.24) 



2/i 



so that if we take fi ~ p, C2n(/^) ~ l/p n+1 . A factor of V 2n at a vertex scales as p 
while each loop contributes a factor of p. Therefore, the leading order contribution to the 
scattering amplitude A-i scales as p" 1 and consists of the sum of bubble diagrams with C 
vertices; contributions to the amplitude scaling as higher powers of p come from perturbative 
insertions of derivative interactions, dressed to all orders by Cq. The first three terms in the 
expansion are 



A 
A x 



-C 



-C 2P 2 



2 ' 



(C 2 p 2 ) 2 M(fi + ip)/4ii 



C 4 p 4 



l + ^iv + ip) 



1 + + 



(2.25) 



where the first two correspond to the Feynman diagrams in Fig. 

Comparing with the expansion of the amplitude eq. (|2.17|) , these expressions relate the 
couplings Cm to the low energy scattering data a, r n : 



S 



A 



A 




+ 



+ 



+ 



FIG. 2. Leading and subleading contributions arising from local operators. 
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CM = T7 



CM 
CM 



Air 
~M 

4ir 



1 


-n 


+ 


1/a 




1 




-n 


+ 


1/a 




1 





2 



1 9 1^1/ 

4 r ° + 2A^ ( ^ + 1/a) 



(2.26) 



M \—fi+ 1/a, 

~ 1/A, these expressions are consistent with the scaling law in 



Note that assuming r n 
eq. (pi). 

This power counting relies entirely on the behavior of C , 2n(A i ) as a function of fi given 
in eq. (|2.24j ). The dependence of C 2n (/x) on \i is determined by the requirement that the 
amplitude be independent of the arbitrary parameter \i. The physical parameters a, r n enter 
as boundary conditions on the RG equations. 

The beta function for each of the couplings C2„ is defined by 



@2n = 11 



dC 



2n 



d/i 



(2.27) 



and they can be computed by requiring that any physical quantity (e.g. the scattering 
amplitude) be independent of \i. In the PDS scheme, the \i dependence of the Cm coefficients 
enters either logarithmically or linearly, associated with simple 1/(D — 4) or 1/(D — 3) poles 
respectively. The functions /3 2n follow straightforwardly from fi-^(l/A) = 0, using the 
expression for A in eq. (|2.23| ). This gives 



0: 



2/i 



Mil 
An 



^ Ci m C2{n- 



in ) 



(2.28) 



m=0 



These /3-functions can also be computed from the one-loop diagrams shown in Fig. 
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ft,: 




XX 



XX 



FIG. 3. Graphs contributing to the ^-functions for C- 



2n 



We examine the RG equations for the first two couplings, Co and C2, in order to explicitly 
show how one recovers the results in eq. ( [2.26|) from solving the renormalization group 
equations. From eq. (|2.28|) it follows that 



Mu 2 

~A t ~ y| 

47T 



• 



4-7T 



CqC- 



2 ■ 



(2.29) 
(2.30) 



Integrating these equations relates the C^n coefficients at two different renormalization scales 
fi and yU - Comparing the theory with A and its derivatives at [A — p — determines the 
initial values C 2ra (0) as in eq. fl2.15p . The solution for C (fi) is 



CM 



CoG^c 



1 + C (no)M(n - /i)/47T 



(2.31) 



With the boundary condition at fi = provided by eq. ( |2.15| ), C (0) = 47ra/M, we arrive 
at the result derived previously for C (/i) in eq. ( |2.26| ): 



° M ~ Tl (-/x+l/a, 



(2.32) 



The RG equation for C2 yields 



C 2 (ji) = C 2 (/i ) 



Co(fio) 



[2.33) 



which when combined with the value of Cq(0) and the boundary condition, C 2 (0) 
Co(0)ar /2, from eq. ( |2.15|) , yields for C^/z) 
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as we found previously in eq. (|2.26|) . 

It is possible to solve the complete, coupled RG equation 



CW — C2mC2( n -m) (2.35) 

for the leading small n behavior of each of the coefficients C2 n ■ The solution, for n > is 

CM = M( \ : , v f i ) +0(^ n ). (2.36) 

M{—fi + l/a) \— im + l/a J 

First note that the scaling property in eq. (|2.24j) is realized: C 2n (/i) oc yU~( n+1 ) for /i ^> \l/a\. 



What is curious is that this leading behavior does not entail a new integration constant for 
each n, but only depends on the two parameters a and r$ encountered when solving for 
Co(/i) and C2(a0; this is due to a quasi-fixed point behavior of the RG equations — the C*2 n 
couplings are being driven primarily by lower dimensional interactions. One can see this 
explicitly in our formula eq. (|2.26|) for C 4 , where the leading 0(/x~ 3 ) part of C 4 depends only 
on r , while the subleading 0(/i~ 2 ) part is proportional to r%. 

This behavior allows us to establish a connection between the present work, and the 
method of introducing an s-channel dibaryon discussed in ref. [§]. The leading /i behavior of 
all of the C2n coefficients is determined by the effective range r^. If one resums this leading 
behavior at the NN vertex one finds (for fi ~ p ^> l/|a|) 

£ c 2n (/V n = S t + o( P 2 ) 

(Stt/MVq) 

E-(-ti+l/a)/Mr ' 1 ' ' 

This looks like an s-channel propagator for a particle at rest of mass [2M + (— fi + l/a)/Mr ], 
and in fact, for fi = 0, corresponds exactly to the dibaryon proposed in || to reproduce the 
scattering due to a short range potential. We see that using the dibaryon is as good as (but 
no better than) carrying out the effective field theory calculation to 0(p l ). The subleading 
corrections can be accounted for by including the subleading part of the C±{ix)p A vertex 
proportional to r±, and which occurs at 0(p 2 ). This dibaryon was recently used with great 
success in the three-body problem 



C. Luke-Manohar velocity scaling 



In ref. |12| a rescaling of fields was shown to eliminate the heavy mass scale M from the 
Lagrangian, and allow a simple power counting in terms of the velocity v of the interacting 
nucleons. One finds that a four fermion operator scales as v , and is therefore irrelevant 
as v — > 0. It is instructive to see why this conclusion is avoided in the theory with large 
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scattering length, where the interaction has a large effect on scattering at the low momentum 

p ~ i/H- 

The rescaled theory is described in terms of dimensionless energy and momentum vari- 
ables, £ = E/Mv 2 and V = p/Mv, corresponding to a rescaling of the spacetime coordinates 
T = tMv 2 , X = xMv. The nucleon fields are rescaled as $ = N / (Mv ) 3 ^ 2 , in order to give a 
canonical energy term. Performing the Luke-Manohar scaling on our toy theory in a naive 
manner yields the action 

S = J d 3 XdT & Ud T + -V 2 x) * + C M 2 v{^^) 2 + ... (2.38) 

This demonstrates the familiar result that a contact interaction (5-function potential) is an 
irrelevant operator in nonrelativistic scattering in D = 4 dimensions, since the interaction 
vanishes as v — > 0. This is certainly true when the underlying theory is perturbative, but 
it fails when it is strongly coupled and there is a large scattering length. In that case we 
should be using the running coupling constants discussed in the previous section, keeping in 
mind that the renormalization scale /x must also be rescaled as ft = fi/Mv, since it is treated 
on the same footing as momenta. Replacing Co in the above expression with the value of 
the running coupling Co(/i) from eq. ( |2.32| ) we find the rescaled interaction 

SXdT Ur , (M/tM/) 2 + ... (2.39) 

-fiMv + l/a y ' v ; 

where ft ~ 1. Now it is evident once again that the limit a — > oo corresponds to tuning Co to 
a UV fixed point. Near that fixed point, the four fermion interaction is a relevant operator 
(it doesn't decrease as v gets smaller), which is why it has a large effect on scattering at 
low velocity. Only for the very low velocities v < 1/Ma does it become irrelevant, at which 
point perturbation theory is once again justified. 



III. THE EFFECTIVE FIELD THEORY EXPANSION FOR REALISTIC NN 

SCATTERING 

We now turn to the problem of interest, realistic low momentum nucleon-nucleon scat- 
tering. The goal of an effective field theory treatment is to provide a systematic expansion 
for computing iViV phase shifts at low momenta, which means that the phenomenology of 
the two nucleon system can be approximated to any desired accuracy. While an effective 
field theory treatment will never in practice compete with the extremely precise phenomeno- 
logical potential model descriptions of the iViV phase shifts, it will provide a framework to 
compute relativistic effects and inelastic processes, as well as to reliably estimate errors at 
any order in the expansion. Furthermore, as it is a perturbative expansion, the effective field 
theory approach will hopefully be useful for many-body systems for which potential models 
are computationally intractable. 

Realistic iViV scattering is more complicated than our toy model for several reasons: the 
inclusion of explicit pion fields, the spin and isospin degrees of freedom, and the need to 
include relativistic corrections in the power counting scheme. Nevertheless, we will show 



that the momentum expansion for this theory closely follows the discussion in §11 B, and 
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that for p < m w the expansion improves as the pion mass gets smaller. This is an important 
result, as it means that there is a region of QCD parameters for which our results hold with 
high accuracy. (We do not, however, take the pion mass so small that am^ is also small.) 

Since we are performing a momentum expansion, it is important to understand the scale 
at which it fails. In the 1-nucleon sector, the expansion is in powers of (m n /A x ) and (p/A x ), 
where A x = Anf = 1.6GeV, (/ = 132 MeV being the pion decay constant). However, we 
find that the chiral expansion for iVW scattering entails the new scale, 

A NN = = 300 MeV • (3- 1 ) 

It is disturbing that i] w = (m^/ Ann) = 0.46, suggesting that our expansion will not converge 
rapidly. However, if one looks at the Yukawa piece of the one pion exchange (OPE) potential 
in the ^So channel, it only binds nucleons for i] n > 1.7, while for the true value f] n = .46, 
the NN phase shift from the OPE Yukawa potential never gets larger than ~ .25 radians. 
Thus there is empirical evidence that the chiral expansion parameter for NN scattering is 
about 30%, and it is reasonable that for momenta p ~ m n the expansion will converge fairly 
quickly. 



A. The chiral expansion 

We need to consider a theory of nucleons interacting with pions and with themselves 
through contact interactions consistent with chiral symmetry. Nucleons are described by 
isodoublet fields N, while the pions are described by the field 

^)=e^, H=f^ (3.2) 

-73/ 

and / is the pion decay constant normalized to be 

/ = 132 MeV . (3.3) 
Under SU(2) L x SU(2) R chiral symmetry the fields transform as 

f _> tft = u£ R \ ? £ = £2 _^ LS/? t ; N^UN , (3.4) 

where L, R are constant SU (2) matrices and U (x) is a pion-dependent SU (2) matrix. The 
vector and axial-vector pion currents are given by 

v, = + fdpt) , a, = i -(td„e - (3.5) 

where the axial current and the chiral covariant derivative = (<9 M + V^) transform 
linearly as 

An -> UAptf , Dp -> UDjtf . (3.6) 
In the 1-nucleon sector, the chiral Lagrangian is then given by 
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£2 f2 

£ = N\iD + D 2 /2M)N + — Tr d^d^E + —uoTr M g (S + S f ) + g A N ] A ■ aN + . . . (3.7) 

8 4 

where M q is the quark mass matrix diag(m u , m^), u has dimensions of mass with m 2 = 
ui{m u + m^), and the ellipses refers to all additional operators. 

When considering the 2-nucleon sector, one must include local four nucleon operators 
in the effective lagrangian. Due to spin and isospin degrees of freedom, the number of four 
fermion operators grows rapidly with the number of derivatives. To the order we will be 
working, we need only consider operators without pion fields and with < 2 derivatives, or 
no derivatives and one insertion of the quark mass matrix. Rather than writing down the 
operators in the chiral Lagrangian, it is simplest to identify their matrix elements between 
particular partial waves. With no derivatives, these contact interactions only effect S'-waves. 
Thus there are two independent Co operators associated with the x 5o and 3 5i channels [|J. 
Denoting the partial waves as 

\a) = \SLJm,E) (3.8) 

we normalize the two operators by defining the Born amplitude for Cc > t ne four fermion 
interactions involving no derivatives, as 

(a\CcM = - (^) E^M* , 7 = {%, 3 Si} . (3.9) 

There are seven independent two derivative contact interactions Q. These operators 
change orbital angular momentum by AL = or AL = 2. Therefore, in analogy to eq. fl3.9| ), 
we can define seven couplings by 



(a\CcM = -[^ l )p 2 



271, L . 

7 = {XM^PoM,^} ■ (3.10) 



As before, p 2 = ME. 

At the same order as the C2 operators are four fermion operators with a single insertion 
of the quark mass matrix M q . If we ignore isospin violation, then M q cx m 2 times a unit 
matrix in isospin space, and so these operators have the same structure as the 0$^ operators. 
We parametrize them as 

(a\£ D2 \P) = - (^) mlY,D?5 ai 5 M , 7 = {X,^} . (3.11) 



7 



Since the ^So and 3 Si scattering lengths are large compared with l/m n , the power count- 
ing in these channels is the same as for the example of §2: for 11 ~ m^, Cq (fx) oc 
C^ 7) (/x) oc l/ii 2 and D£\n) oc 1/li 2 for 7 = {%, 3 5i}. A single exchange of a potential pion 
contributes to the amplitude a factor of 

2/ q + Tn* 
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FIG. 4. Leading contributions to the [3-functions for Cq , including pions 



which scales as 0(1) in the chiral expansion. Thus one pion exchange occurs at the same 
order as an insertion of C2P 2 or Divri^ in the and 3 5i channels, and pion exchange may be 
treated perturbatively. A new feature of the theory with pions is that this scaling behavior 
breaks down not only at low momentum, p ~ l/\a\, but also at sufficiently high momentum. 
This can be seen by examining the RG equations for C^. The exact beta function for the 
Cq 1 ^ coefficients comes from the graphs of Fig. [|, and one finds 




(3.13) 



^ = ^ = ^ 1 (cfl))2 + 2cfl)| | +9 (||y|. (3 , 4) 



Solving the RG equation for the singlet channel eq. ( |3.13|) with the boundary condition 
Cq S °\o) = 47rai/M, where a\ is the ^So scattering length, we find 

gf-'W = ~ L + A) . (3.15) 

where Ajvjv is given in eq. Since Ajvtv 3> for /i ^> 1 / 1 ax | we have 

and so the power counting changes for /i > Ann; in fact the UV fixed point toward which 
Cq is driven largely cancels the 5-function component of the single pion exchange in the 

channel. Similarly, power counting is found to change at ji ~ K NN in the channel 
as well. As a result, the power counting developed in §2, with one pion exchange treated 
as 0(1), works only up to p ~ Ann- At that point one pion exchange and a Cq insertion 
become equally important. The power counting that we follow in this paper is therefore 
expected to fail completely at momenta on the order of Ann- 

Up to now we have only discussed the scaling of the O2 interactions in the ^So and 3 S\ 
channels. The l//i 2 scaling of these couplings is due to multiplicative renormalization by 
Cq. Such renormalization does not occur for O2 in any of the P-wave channels, and so 
these couplings will be 0(p°). Thus the leading contribution to the P-wave (and all higher 
partial wave) amplitudes occurs at 0(1) and is simply a single pion exchange (the Born 
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FIG. 5. Graphs contributing to the leading amplitude A-\. 



approximation OPE). It is well-known that this accounts very well for the observed phase 
shifts. The subleading contribution to the P-wave amplitude is due to the exchange of two 
potential pions at OQo 1 ), while the C2P 2 operators enter at 0(p 2 ), along with the two- loop 
graph with the exchange of three pions. 

The transition operator coefficient, C 2 Sl ~ is different from the others as it is multi- 
plicatively renormalized by Cq , but its beta function is half as big as that for C 2 . For 
small /1, C 2 Sl Dl \fJ>) ~ and it enters the calculation at 0(p 1 ). An important aspect of 
the OQo 1 ) calculation of the ^So and 3 5i scattering amplitudes is the absence of new undeter- 
mined coefficients C[ S °^ and C[ Sl \ The renormalization group scaling determines the lead- 
ing behavior of these couplings. For small /i (but \a\fi 3> 1), c!f\fi) = — (Mfi/4ir)C 2 (a*) 2 - 

It is now possible to perform a systematic chiral expansion for iViV scattering amplitudes 
in the situation where m n is small, (but |a|m T is not). Denoting our expansion parameter 
by Q, representing equally p or m T , and taking \i ~ Q the leading amplitude is 0(Q~ 1 ) and 
is simply the bubble sum of interactions, shown in Fig. [|. The subleading contribution 
is 0(Q°), and involves an insertion of either a single pion exchange, a C 2 operator, or a D 2 
operator, dressed to all orders by the C interaction, as shown in Fig. |. At 0(Q 1 ) one finds 
both two pion exchange (involving both the two pion vertex as well as iterated one OPE), 
four derivative operators, and so forth. It may be tempting to sum up the OPE potential 
exactly as is typically done (e.g., ref. 0]), but it is not consistent to do so while leaving out 
the arbitrarily higher derivative contact interactions that are required to renormalize the 
pion ladders. By neglecting these higher derivative terms one is making a model for the 
short distance physics, and the result one gets is not any more accurate for having included 
the multiple pion exchange. 



B. The channel to subleading order 

In the isotriplet channel the nucleons are exclusively inanS' = L = J = state. The 

/let \ 

leading 0(Q^ 1 ) contribution to scattering in this channel is from the bubble chain of Cq 
operators, as shown in Fig. || 



C, 







A-i = ^f- . (3.17) 

l + cf 0) g(// + ip) 

Since we are computing the amplitude between plane wave states, we omit from A the 
normalization factor (Mp/2ir) appearing in eqs. (|3.9| - |3.1ll) . There is only one unknown 



parameter that needs to be fit to data and for simplicity we will require that the amplitude 
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at this order reproduce the (np) scattering length a\ = —23.714 ± 0.013 fm. This 
determines CQ S °\m n ) = —3.51 fm 2 , where we have chosen to renormalize at /i = m % and 
the behavior of the phase shift over a range of momenta resulting from this fit is shown in 
Fig. ft 

At next order Q° there are contributions from insertions of higher dimension local oper- 
ators and also from potential pion exchange. The 0(Q°) contribution to the amplitude is 
written as a sum of the five terms, Aq = Aq + Aq 11 ^ + Aq 1 ^ + A + A . The local oper- 
ators at this order involve either two spatial derivatives, C% , or one insertion of the light 
quark mass matrix, D\ . Expressions for the graphs shown in Fig. ^ which we presented 
in ref. 01 are 



A^ 



A^ 



A HI) 



A [IV) 
A (v) 



1 2 




\2f 2 J \ 4p 2 
P V 




—D 



A-i 



tan" 



i tan 



2 ^ + ilnfl + ^ 

2 V ml i 

__ - l n ( m * + 4p2 



i l 2P 



ni 



+ 1 



(3.18) 



The two-loop diagram in Aq V ^ involving the exchange of a potential pion between two 
contact terms is divergent in both three and four dimensions. In the PDS scheme we 
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subtract the poles in three and four dimensions leaving this graph logarithmically as well 
as power-law dependent on the renormalization point fx. As the coefficient of the four- 
dimensional divergence is proportional to the mass of the pion squared, the required isospin 
conserving counterterm with coefficient D 2 (fx) depends on the sum of the light quark 
masses, m u + m d and gives rise to Aq"\ In addition, for convenience we have absorbed 
some of the /^-independent terms from Aq V ^ into the definition of D 2 So \fx). At this order 
there are three unknown counterterms that need to be fit to data, CQ So \fx),C 2 So \fx) and 
D2 S °\n)- As the amplitude can be written as a function of Cq So \ix) + m^D^^dx) , the 
dependence of observables upon C ( Q So) (fx) and D ( 2 1S "\fx) individually is an artifact of the 
perturbative expansion, and is indicative of the size of higher order effects. Conventionally, 
the scattering data in the NN sector is presented in terms of phase shifts. In this channel, 
the phase shift is simply related to the amplitude A by 



S = e 2iS = l + i^-A , (3.19) 



5= ^ln(l+t^A) . (3.20) 



leading to 



2i V 2tt 

Expanding both sides to a given order in Q with 5 = + 5^ + . . gives 

*>-I h ( 1 + #u), ^(iriid- 

Here superscripts denote the order in Q. Our expression for the amplitude A gives an S- 
matrix that is unitary up to the order we have computed, i.e. if the amplitude is computed 
up to 0(Q k ), then S*S = 1 + 0(Q k+2 ). 

It is convenient to choose fx = m 7r . Expressions for the scattering length and effective 
range are determined from the expansion, pcotS = ip + 4it/MA-i — 4iiAo/MA 2 _ 1 , which 
yields to the order we are working, 

«i V MC^Hm,)) M(Ci s "(m,)f 

and 

r = 8nC^\ m7T ) g\M_ ( 16vr 96vr 2 \ 

° M(Cf o) (m.)) 2 l^P\ cS So) (m«)m«M (C™ \m^Mf ) ' 

The choice 

C (m n ) = -3.63 fm 2 , D 2 (m w ) = , C 2 (m n ) = 2.92 fm 4 , (3.24) 
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is consistent with the experimental values of the scattering length and effective range []. 
About 43% of the effective range is due to C^m^)- The phase shift resulting from these 
parameters is shown in Fig. [7]. 

Alternatively, we consider fitting the phase shift from the Nijmegen partial wave analysis 
|32|1 over the momentum range p < 200 MeV, treating Cq So \ D^ ^ and C\ ^ as free 
parameters. The results are 

c£ So \m n ) = -3.34 fm 2 , D^ S °\m^) = -0.42 fm 4 , C^ S(, \m^) = 3.24 fm 4 , 

(3.25) 

which give the phase shift plotted in Fig. 0. As is apparent from Fig. |7|, the agreement 
of the phase shift with data is excellent at quite large values of p. Furthermore, the cou- 
pling Cq (TOtt) is close to its leading order value (in the limit of large scattering length), 
— (Air/Mm^) = —3.7 fm 2 , and C2 (ra^) is also near its expected size, suggesting that our 
expansion is valid in this channel. However, for p > 100 MeV the magnitude of the ratio 
Ao/A-i is greater than ~ 0.5 and it is difficult to justify the approximations we have made, 
e.g. neglecting terms suppressed by (Aq/A-i) 2 . The difference between the two fitting 
procedures is attributable to effects higher order in our expansion. 

C. The 3 Si — 3 D\ channels to subleading order 

The analysis of the isosinglet 3 Si — 3 Di channel is richer than the 1 5o channel since 
there are two different orbital angular momentum states involved. The power counting in 
the 3 5i channel is the same as for the ^So channel. However, we must also consider how the 
coefficients of the operators contributing to scattering in the 3 D\ channel and the coefficients 
of the operators that give rise to 3 «Si — 3 D\ mixing behave under renormalization group scaling 
and at what order in the Q expansion they contribute to observables. Firstly, operators 
between two 3 D\ states are not renormalized by the leading operators, which project out 
only 3 5i states. Further, they involve a total of four spatial derivatives, two on the incoming 
nucleons, and two on the out-going nucleons. Therefore, such operators contribute at 0(Q A ), 
and can be safely neglected. Consequently, amplitudes for scattering from an 3 Di state 
into an 3 Di state are dominated by single potential pion exchange which contributes at 
O(Q ). Single pion exchange describes well this partial wave in the momentum range we are 
considering. Secondly, operators connecting 3 D\ and 3 S\ states involve 2 spatial derivatives 
(acting on the 3 D\ state) and are renormalized by the leading operators, but only on the 

(3S — 3 D) 

L = "side" of the operator. Therefore the coefficient of this operator, G\ 1 ~ l//x, 
and contributes at OiQ 1 ). Hence it can also be neglected at the order we are working. Thus, 
mixing between 3 Di and 3 <Si states is dominated by single potential pion exchange dressed 



At this order there is ambiguity in these values since to the order we are working only the linear 
combination, Co(?n 7r ) + m^D2(m. K ), is determined. At higher order the Co and D2 operators are 
distinguished by a contribution to an NN^ Ntttt vertex proportional to D2. 
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p (MeV) 

FIG. 7. The phase shift 5 for the 1 So channel. The dot-dashed curve is the one parameter fit at 
0(Q^ 1 ), that reproduces the scattering length. The dotted and dashed curves are the fits at O(Q ) in 
eqs. ({3.24 , 3.23j ) respectively. The dashed curve corresponds to fitting 5 between < p < 200 MeV, 
while the dotted curve corresponds to fitting the scattering length and effective range. The solid line 
shows the results of the Nijmegen partial wave analysis. 



by a bubble chain of Cg operators and a parameter free prediction for this mixing as a 
function of momentum exists at O(Q ). 

We denote the amplitude at 0{Q n ) by A n [LL'}, where L and L' are the initial and final 
orbital angular momenta. As in the ^So channel, we omit from A the normalization factor 
(Mp/27i) appearing in eqs. (|3.9| - |3.11| ) since we are computing the amplitude between plane 
wave states. At leading 0{Q~ 1 ) in the expansion there is a contribution only to the 3 S± 
partial wave: 



r ( 3 Si) 



A 



-1[02] 



A-\ [20] 



A- 



1[22] 



(3.26) 



At O(Q ) there are contributions from graphs of the same form as in the amplitude for 
scattering, shown in Fig. [!]. Using the same identification of graphs as in the ^So channel, 
and the similar notation, Ao[l,l'] = A^) L L n + ...,we find that 



A 



'A-i\, 



a 



( 3 Si) 



A (I) 

,Al 0[02] 



A {1) 

,Al 0[20] 



A {1) 

,Al 0[22] 



(3.27) 



arising from a single insertion of the local operator involving two spatial derivatives. Single 
potential pion exchange in the t-channel gives 



A 



in) 



2/2 



1 



fry 



20 



A in) 



A 



in) 

0[22l 



(II) . 
0[20] 

9a™1 
4f 2 p 2 



9a 



V2f 



[Qo{z) + Q 2 {z)-2Q l {z)\ , 



6p 2 



Q2(z) + -^(Q 1 (z)-Q 3 (z)) 



(3.28) 



where z — 1 + ml / \2p 2 ) and Qk(z) denotes the fc-th order irregular Legendre function, 

Q ° = ^ l0g (^l) ' Q^ z Qo~i, (n + l)Q n+1 = (2n + l)zQ n -nQ n . 1 . (3.29) 

The contribution from single pion exchange across the end of a bubble chain of operators 
with coefficient Cq Si ^ is 



(in) _9a( m n MA-i[, 
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^ f M„4_l[ 00 ] \ ^ 2 



47T 
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(/v) _ 9a ( m^MA-iioo] 



•A) [00] 
*^0[02] 
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,/H -0[20] 
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fi + ip 
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is from pion exchange between two chains of operators with coefficients C . Finally, a 
single insertion of the quark mass matrix leads to 



(3.30) 

■)■ 

(3.31) 



A)[00] ~~ u 2 m ir 



U 



A (v) 

,Al 0[02] 



A 



(V) 
0[20] 



A iv) 

,Al 0[22] 







(3.32) 



Again part of the subtraction point independent contribution to -4 [oo] nas b een absorbed 
into A { ^ y 

The S-matrix in this channel is usually expressed in terms of two phase shifts, 5o and 62, 
and a mixing angle £1, 



5 = 1 + i P -^A 



e 2i5 ° cos 2ei i e ^o+s 2 ) sin 2£l 
ie i(s +s 2 ) sin 2ei e 2i<52 cos 2e 1 



(3.33) 



and like the channel we will expand S order by order in Q. 

As 3 <Si — 3 Di mixing has vanishing contribution at order 0{Q~ r ) the mixing parameter S\ 
starts at 0{Q 1 ) ) the same holds true for 5 2 (the explicit factor of p in the relation between 
S and A increases the order by 1). Writing each of the parameters as an expansion in Q, 
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So = 5$ 0) + + - , 6 2 = 4 0) + + - , ex = 4 0) + + ... , (3.34) 



it follows that 
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log 



1 + 1 17" 4 - 



1[00] 



x(l) 
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0[00] 



47T 



1 + A. ir 



(3.35) 



4 0) = , e« = , A[02] , (3.36) 

4-^/TT^Xx^' 

4 0) = , = ^ Am . (3.37) 

Working to subleading 0(Q°) there are three parameters describing the 3 Si phase shift. 
Again it is convenient to choose the subtraction point equal to m^. However, as we dis- 
cussed previously observables do not depend upon CQ Sl \m^) and D^^im^) indepen- 
dently. Therefore we can set D^m^) = and fit Cq and C 2 to the low energy observables 
taken to be the scattering length a 3 = 5.423 ± 0.005 fm and the deuteron binding energy 
E D = 2.224644 ± 0.000034 MeV. The result of this fit is 

c£ Sl) (m w ) = -5.03 fm 2 , D ( 2 Sl \m n ) = , c£ Sl) {m n ) = 4.96 fm 4 . (3.38) 

Alternatively fitting the parameters Cq Sl \m n ), C 2 Sl \m 7r ) and D^^m^) to the phase 
shift So over the momentum range p < 200 MeV yields 

c£ Sl \m n ) = -5.51 fm 2 , D { 2 Sl) {m n ) = 1.32 fm 4 , cf 5l) (m^) = 9.91 fm 4 . (3.39) 



Fig. ^ shows the phase shift compared to the Nijmegen partial wave analysis |32[ for this 
latter set of parameters. The other two quantities, S\ and S 2 , receive no leading order 
contributions and both begin at O(Q ). There are no free parameters at this order in either 
Si or 5 2 once C has been determined from Sq. The predictions for e\ and S 2 from the fit 
eq. (|3.39|) and a comparison to the Nijmegen partial wave analysis can be found in Fig. |8|. 



D. Higher partial waves 



The analysis of the previous section demonstrates how the power counting impacts the 
3 D\ channel and this discussion generalizes to other partial waves. A local operator that 
connects an angular momentum L state with an angular momentum V state involves at 
least L + L' spatial derivatives. The case of S'-wave to S'-wave scattering has been described 
in the previous sections. If either L or V but not both correspond to an S'-wave then the 
operator enters at 0(Q L+L However, if neither L nor V is equal to zero the operator 
contributes at 0(Q L+L ') for L, L' odd and O^Q 1 ^'- 1 ) for L, L' even. The contribution of 
pions is at O(Q ), and is therefore the leading contribution to all non S'-wave to S-wave 
scattering amplitudes. This contribution has been presented in the literature (e.g. ||33|| ) . 
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FIG. 8. The phase shifts So, 82 and mixing parameter e\ for the 3 S\ — 3 D\ channel. The solid 
line denotes the results of the Nijmegen partial wave analysis. The dot-dashed curve is the fit at 
0(Q~ 1 ) for 5q, while 62 = £\ = at this order. The dashed curves are the results of the 0(Q°) fit 
of 5q to the partial wave analysis over the momentum range p < 200 MeV, as given in eq. I[3. 3(\) . 



E. Radiation pions and operator mixing 



We have seen that graphs involving potential pion exchange occur at O(Q ) in both 
the ^So and 3 S\ — 3 Di channels. Such contributions arise from kinematic regions where the 
intermediate nucleons are near their mass-shell, with the pion exchanged in the t-channel far 
from its pole. Contributions arising from radiation pions (the pion pole) exchanged in the 
s,t and u channels arise at 0(Q l ) (for discussion of radiation exchanges in nonrelativistic 

An interesting feature of virtual radiation pions in the s-channel 



gauge theories see p6H31[] 
is that they can cause mixing between four-nucleon operators in different spin channels. This 
is because the two nucleons in a virtual intermediate NNir state can rescatter while in a 
different isospin and angular momentum state than the physical incoming nucleon pair. 

As an example, consider the graph shown in Fig. |9| with an insertion of the operator 
(N^a a N) 2 , denoted by IV Explicitly, 



Ti = i^-N ] T a a i a k NN ] T a a k a i N 

4-D , d D q 
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+ ie 



. (3.40) 



The q° integral is performed by forming a contour enclosing the one pole in the upper 
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FIG. 9. A graph involving a radiation pion that gives rise to mixing between the So — 3 S*i 
operators. This graph with an insertion of (N^a a N) 2 at the vertex gives the expression F\ in 
eqs. KMKM. 



half of the complex plane provided by the pion propagator. Using the equations of motion 
E = p 2 /M + 0(p 4 /M 3 ) and neglecting the 0(p 4 /M 3 ) relativistic correction we realize that 
the weight of the integral in the low energy theory will be for momentum near the pion mass 
and perform an expansion in 1/M, giving 



r\ = ^N^T a <T i a k NN^T a a k a i N (^j "J 
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+ 0{1/M) . (3.41) 



Evaluating the integral yields 
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Including all the irreducible graphs and wavefunction renormalization it is straightfor- 
ward to find the leading radiation pion contribution to the (3 functions for D\ and D 2 

are 
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(3.43) 



These contributions give rise to mixing between the S"-wave spin-singlet and spin-triplet 
operators. 



F. Relativistic effects 

A further contribution starting at 0{Q X ) arises from relativistic corrections to the energy- 
momentum relation. A detailed discussion of such effects in dimensionally regulated non- 
relativistic gauge theories can be found in |29[ and the situation is similar for nucleon 
interactions. Neglecting pion fields, the lagrange density in the single nucleon sector is 

c = N< ( !S ° + m) N + Mi N>v ' N + - ■ (3 - 44) 

where the spatial gradient operator brings down factors of p. It is understood that the 
N Jf 'V 4: N operator is inserted perturbatively into graphs, e.g. Fig. [TO, and that the lowest 
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FIG. 10. Higher dimension operators arising from relativistic corrections to the en- 
ergy-momentum relation, denoted by the crossed- circle, are inserted perturbatively. 

order equations of motion are modified to E/2 = p 2 /2M + p 4 /8M 3 + . . .. A single insertion 
of N^V^N gives rise to a 0{Q l ) contribution; however, it is suppressed by factors of the 
nucleon mass and not A NN . Consequently, its effect is expected to be small compared to 
other corrections at this order. 

IV. OVERVIEW 

In a previous letter we presented a new power counting scheme to describe NN scat- 
tering processes that does not suffer from the inconsistencies of Weinberg's scheme [||. In 
order to achieve consistent power counting a new subtraction scheme is used for dimen- 
sionally regulated integrals, designed to keep track of power law divergences in theories in 
which there are delicate cancellations at short distance; for NN scattering, these cancella- 
tions manifest themselves as a large scattering length. The renormalization group provides 
a powerful tool in the analysis of such theories and allows one to identify the order at which 
a given operator will contribute. 

In the present paper we have elaborated on our expansion, with particular emphasis on 
the utility of the renormalization group. We have also presented detailed analytic compu- 
tations for NN scattering in the ^So and 3 S\ — 3 Di channels, complete to the subleading 
0(Q°) order, and have shown that the agreement with experiment at low energy is quite 
good. A particular success was the calculation of the 3 S\ — 3 Dx mixing parameter E\ shown 
in Fig. H, with no free parameters. We then discussed the power counting for higher partial 
waves, and showed that they are dominated by pion exchange in the Born approximation. 
Finally, we showed how to deal with relativistic effects and virtual radiative pions, which 
are features that arise at order OiQ 1 ). 

The techniques presented here should be applicable to a number of low energy processes, 
such as radiative capture, and electromagnetic moments of the deuteron. Theoretical chal- 
lenges include extending the validity of the expansion above the scale A^at, and applying 
the technique to systems with three or more nucleons. 
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